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ABSTRACT 

A new classification scheme is introduced for comet-likc bodies in the Solar system. 
It covers the traditional comets as well as the Centaurs and Edgeworth-Kuiper belt 
objects. At low inclinations, close encounters with planets often result in near-constant 
perihelion or aphelion distances, or in perihelion-aphelion interchanges, so the minor 
bodies can be labelled according to the planets predominantly controlling them at 
perihelion and aphelion. For example, a JN object has a perihelion under the con- 
trol of Jupiter and aphelion under the control of Neptune, and so on. This provides 
20 dynamically distinct categories of outer Solar system objects in the Jovian and 
trans- Jovian regions. The Tisserand parameter with respect to the planet controlling 
perihelion is also often roughly constant under orbital evolution. So, each category can 
be further sub-divided according to the Tisserand parameter. 

The dynamical evolution of comets, however, is dominated not by the planets 
nearest at perihelion or aphelion, but by the more massive Jupiter. The comets are 
separated into four categories - Encke-type, short-period, intermediate and long-period 
- according to aphelion distance. The Tisserand parameter categories now roughly cor- 
respond to the well-known Jupiter-family comets, transition-types and Halley-types. 
In this way, the nomenclature for the Centaurs and Edgeworth-Kuiper belt objects 
is based on, and consistent with, that for comets. Given the perihelion and aphelion 
distances together with the Tisserand parameter, our classification scheme provides a 
description for any comet-like body in the Solar system. The usefulness of the scheme is 
illustrated with examples drawn from numerical simulations and from the present-day 
Solar system. 

Key v^rords: minor planets, asteroids - comets: general - planets and satellites: 
general - celestial mechanics - Kuiper belt - Solar system: general. 



1 INTRODUCTION 



There is a need to improve the taxonomy for comet-like 
bodies in the Solar system. First, as a result of recent dis- 
coveries, the populations of comet-like bodies in the Solar 
system are known to be much more extensive than pre- 
viously thought. The last decade has seen the discovery 
of more than 600 trans-Neptunian objects, beginning with 
(15760) 1992 QBi (Jewitt & Luu 1993), as well as the iden- 
tification of ~100 Centaurs (e.g. Scotti 1992) to supple- 
ment the early serendipitous finding of (2060) Chiron (Kowal 
1979). Moreover, many unusual objects have now been found 
circulating essentially within the known planetary system, 
for example (5335) Damocles (e.g. Asher et al. 1994), and 
the 'asteroids' (4015) Wilson-Harrington = 107P/Wilson- 
Harrington, (7968) Elst-Pizarro = 133P/Elst-Pizarro and 
(2060) Chiron = 95P/Chiron. These objects blur the tradi- 



tional clear distinction between comets and asteroids, either 
on physical or dynamical grounds. 

Secondly, the spatial distribution of objects in orbits 
largely beyond Neptune, but with aphelion distances Q 
much less than the conventional inner edge of the dynami- 
cally active Oort cloud (Q<c20 000au, say), is believed to 
comprise a flattened disc or belt-like structure. These objects 
have orbits that in some cases extend hundreds of astronom- 
ical units (au) from the Sun, and although resembling aster- 
oids in appearance are widely thought to be cometary in 
composition. This trans-Neptunian zone includes the region 
referred to as the Kuiper or Edgeworth-Kuiper belt (EKB), 
and probably merges in the form of an extended trans- 
Neptunian disc into the inner core of the Oort cometary 
cloud. Trans-Neptunian objects are a possibly significant 
source of low-inclination 'Jupiter-family' comets (Quinn, 
Tremaine & Duncan 1990) , again muddying the traditional 
distinction between comets and asteroids. 
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In fact, the Centaur population appears to represent 
an important link between trans-Neptunian objects and 
Jupiter-family comets (Dones, Levison & Duncan 1996, 
Stern & Campins 1996, Levison & Duncan 1997). Centaur 
orbits arc typically planet-crossing and have relatively short 
dynamical lifetimes (~10^yr). Chiron, which is one of a 
number of exceptionally large minor bodies with perihelia 
close to or within the orbit of Saturn, exhibits cometary 
activity (e.g. Luu & Jewitt 1990) and even has a periodic 
comet designation: 95P/Chiron. At least one recently dis- 
covered comet, namely C/2001 T4 (NEAT), has a very sim- 
ilar Chiron-like Centaur orbit. Taken together, the evidence 
suggests a picture in which comets and distant minor plan- 
ets arc dynamically reprocessed from one population to an- 
other, for example from the Edgeworth-Kuiper belt, through 
Centaurs to Jupiter-family comets, and so on. Similarly, it 
is possible for objects to be scattered out of the main as- 
teroid belt and evolve into similar areas. So, it is possible 
that some Centaurs may be rocky or astcroidal, others may 
be icy or cometary. This is analagous to the examples given 
by Fernandez, Gallardo & Brunini (2002) of mixing between 
the Jupiter-family and near-Earth asteroid populations, and 
again highlights both the difficulty of separating physically 
distinct populations of objects using purely dynamical char- 
acteristics, and the need for a unified dynamical classifica- 
tion scheme capable of describing the full range of observed 
orbital types. 

In parallel with these observational advances, the march 
of computer power now ensures that the orbits of most So- 
lar system bodies can be routinely integrated for millions of 
years. Given the wealth of simulation data and the diver- 
sity of new discoveries, the taxonomy of Solar system ob- 
jects assumes great importance. For a classification scheme 
to be useful, it should allow us to place objects with similar 
physical or dynamical characteristics in sets, and to examine 
the bulk statistics of the sets in detail. However, in making 
sense of the complicated evolutionary pathways followed by 
objects such as short-period comets and Centaurs, the his- 
torical classification scheme seems obsolete. A new taxon- 
omy should clarify the principal dynamical paths followed 
by different classes of object through interplanetary space, 
and facilitate the definition of dynamical lifetimes and the 
flux and transfer probabilities of objects from one dynamical 
class to another. 

Traditionally, comets have been classified as short- 
period or long-period according to whether their periods of 
revolution P are less than or greater than 200 years. This 
unphysical division (Levison 1996) has frequently been jus- 
tified in terms of the time over which comets have been ob- 
served with scientific precision (e.g. Weissman 2001). But, 
the canonical upper limit for 'periodic' comets has now been 
spectacularly broken by the return of 153P/Ikeya- Zhang 
with an orbital period in excess of 350 years, and also by 
recent discoveries such as asteroid 2002 RP 120 with a period 
of about 420 yr. 

Centaurs are conventionally defined as asteroids or 
cometary nuclei circulating largely between the orbits of 
Jupiter and Neptune and usually crossing the orbit of at 
least one giant planet (e.g. Jedicke & Herron 1997, Larsen 
et al. 2001). Jedicke & Herron's definition restricts Centaur 
semi-major axes to less than that of Neptune (i.e. 0;$ 30 au), 
and therefore allows objects defined as Centaurs to pass 



significantly beyond Neptune's orbit into the Edgeworth- 
Kuiper belt. A straightforward extension of this definition 
would allow Centaurs to be regarded simply as a continu- 
ation of the 'scattered disc objects', perhaps even encom- 
passing those resonant objects that are not protected in the 
long term from close approaches to Neptune. The scattered 
disc was the name introduced by Duncan & Levison (1997) 
following the discovery of 1996 TLee (Luu et al. 1997) to 
describe a population of relatively high-inclination, high- 
eccentricity trans-Ncptunian objects. The idea that Cen- 
taurs and the scattered disc may be drawn from the same 
underlying population has been advocated recently both 
by Marsden (1999a) and by Emel'yanenko, Asher & Bailey 
(2003). 

This paper introduces a new classification scheme for 
comet-like bodies of the Solar system. There is much ev- 
idence from previous work (e.g. Kazimirchak-Polonskaya 
1972; Kresak 1972, 1980, 1983; Everhart 1977; Rickman & 
Frocschle 1988; Manara & Valsecchi 1991; Dones et al. 1996; 
Levison 1996) that the evolution of comets and Centaurs of- 
ten takes place under the control of one or another major 
planet, usually with near-constant aphelion or near-constant 
perihelion, or via aphelion-perihelion interchanges. This mo- 
tivates us to consider the planets controlling the perihelion 
q and aphelion Q as fundamental to the classification of any 
comet-like bodies. There is also evidence that low-inclination 
objects tend to evolve inwards, whereas the reverse is true 
for high-inclination objects. This suggests that we use the 
inclination - or better still, the Tisserand parameter Tp with 
respect to the controlling planet - as a third criterion. The 
Tisserand parameter is an approximate constant of the mo- 
tion, at least when the orbital perturbations are dominated 
by a single planet. Given the set of values (q,Q,Tp), our 
new taxonomy provides an instantaneous classification for 
any comet-like body in the Solar system. 



2 AN OUTLINE OF THE CLASSIFICATION 

SCHEME 

2.1 The zone of control 

The Hill radius is defined as 



where Mp is the mass of the planet, Mq is the mass of the 
Sun and a is the semi-major axis of the orbit of the planet 
about the Sun (e.g. Murray & Dermott 1999). It corresponds 
to the position of the Lagrange points in the restricted three- 
body problem (e.g. Arnold et al. 1987) and so marks the 
largest distance at which the planet may possess a moon. 
A Solar system object can be classified according to the 
planets which have the greatest effect at perihelion and at 
aphelion. For typical eccentricities, a perihelion lying 4 or 5 
Hill radii beyond an inner planet is often sufficiently distant 
for the dynamics to be controlled by the outer planet. So, 
we conclude that for a planet to control an object, then the 
perihelion or aphelion must lie closer than a distance of «3 
Hill radii to the planet's orbit (see e.g., Charnoz, Thebault 
& Brahic 2001). 
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Tisserand 
Class 



Tisserand 
Parameter 



Class I 
Class II 
Class III 
Class IV 



Tp < 2.0 
2.0 < Tp < 2.5 
2.5 < Tp < 2.8 
2-8 < Tp 



Table 1. Definition of Tisserand parameter classification scheme. 
The boundary near Tp ~ 2.8 marks the limit above which it is im- 
possible for an object to be directly ejected in a single encounter. 



In our classification scheme, each planet's zone of con- 
trol extends to 3k Hill radii, with the parameter fc ~ 1. In 
this way, the Solar system is banded into concentric zones of 
control of first Jupiter, then Saturn, Uranus and Neptune. 
Over the next few million years the mean semi-major axes 
and max;imum aphelia of the four major planets Jupiter, Sat- 
urn, Uranus and Neptune are (5.2, 5.5), (9.6, 10.4), (19.2, 
20.8) and (30.1, 31.0) au respectively (Applegate et al. 1986). 
Three times their Hill radii are 1.1, 1.3, 1.4 and 2.3 au respec- 
tively. These distances provide a rough basis for estimating 
the respective boundaries between the zones of control. For 
example, for Jupiter to control the object at perihelion, the 
perihelion must lie between 4 and 6.6 au approximately. If 
the aphelion also lies in the range, then we classify the mi- 
nor body as a J object. For Saturn to control the object at 
aphelion, then 6.6 < Q < 12.0 when the object is denoted by 
JS. We will give detailed examples shortly, but the general 
idea should now be obvious to the reader. 



2.2 The Tisserand parameter 

The Tisserand parameter is an approximation to the Jacobi 
constant, which is an exact integral of motion in the circular 
restricted three-body problem. It is defined as (e.g. Murray 
& Dermott 1999) 



ap „ . / a(l — t 

h 2 cos i\ — 

a V Hp 



(2) 



where ap is the planet's semi-major axis. 

Ignoring the positive orbital eccentricity of the planet, 
objects with Tp greater than 3 are in principle confined to re- 
gions wholly interior or wholly exterior to the planet's orbit. 
There is a long history of usage of the Tisserand parameter 
with respect to Jupiter Tj to divide the cometary families 
(e.g., Kresak 1972, 1980, 1982, 1983, 1985; Vaghi 1973 a,b; 
Carusi & Valsecchi 1987; Levison & Duncan 1994; Levison 
1996). Fig. Q shows all short-period comets in the plane of 
semi-major axis a and eccentricity e. The comets have been 
colour-coded according to their Tisserand parameter with 
respect to Jupiter. Motivated by these plots, we propose a 
four-fold division according to Tisserand parameter as indi- 
cated in Table 

In our classification scheme, we use the Tisserand pa- 
rameters with respect to the major planets (rj,Ts,ru and 
Tn with obvious notation) to classify objects whose peri- 
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Figure 1. Plot of eccentricity versus semi-major axis for short- 
period comets in Marsden & Williams (1999). All objects are 
colour-coded according to their Tisserand parameter classes with 
respect to Jupiter. Most of the comets lie in the SP category, 
leftward of the q = 4.0 au line. To the right of that line, the 
boundaries of the J, JS and JU categories are marked. 



helion lies within the zone of the planet's control. So, for 
example, an SUiv object is a minor body whose perihelion 
is controlled by Saturn and aphelion by Uranus and whose 
Tisserand parameter with respect to Saturn is larger than 
approximately 2.8. Here and henceforth, we always denote 
the Tisserand parameter class as a subscript. 



3 COMETS 

3.1 Historical introduction 

Objects classified as comets display the widest variety of 
dynamical behaviour in the Solar system. Unfortunately, 
there is little or no agreement as to the precise definition of 
different cometary types. Historically, the major cometary 
division has been between those classified as 'long-period' 
(P > 200 yr) and 'periodic' or 'short-period' (P < 200 yr) 
comets (e.g. Duncan, Quinn & Tremaine 1988, Nakamura & 
Kurahashi 1998, Levison 1996, Crovisier 2001). The major- 
ity of the short-period group belong to the so-called Jupiter- 
family, most of which have very short periods indeed (a 
mean value close to 8 years). For this reason, the periodic 
comets are usually divided simply on the basis of orbital pe- 
riod into two classes: Jupiter-family comets with P<20yr 
(e.g. Fernandez 1985, 1994), and Halley-type comets with 
20<P<200yr (e.g. Carusi et al. 1987b). 

Since 1999, however, one-apparition comets with peri- 
ods in the range 30 < P < 200 yr have been given a standard 
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cometary 'C/' designation, rather than the usual periodic 
'P/' identifier (Marsden 1999b, Green 2000). The official 
range for short-period comets therefore now extends to only 
30 yr instead of the former 200-year limit, and comets with 
30 < P < 200 yr are now described as 'intermediate period' 
(Marsden & Williams 1999). The principal disadvantages 
of this change are that one-apparition intermediate-period 
comets cannot be easily identified as periodic simply from 
their designation, while the phrase 'intermediate period' in 
fact has a rather wide range of previous connotations, hav- 
ing been used to describe objects with periods both in the 
Halley-type 20 < P < 200 yr range (e.g. Fernandez 1980, 
Fernandez & Gallardo 1994) and with much longer periods 
up to ~10^yr (e.g. Everhart 1974). 

In fact, the situation is more confused still. Whip- 
ple (1978) appears to have chosen P = 25 yr to separate 
his 'Class V short-period comets (essentially the classical 
Jupiter-family) from his Class IV intermediate-period group 
(25<P<1000yr); while Everhart (e.g. 1972, 1973) adopted 
13 years instead of 20 years as a convenient boundary be- 
tween short-period and intermediate-period comets, the lat- 
ter being defined to have periods in the range 13 < P < 
1000 yr. Others (e.g. Rickman & Froeschle 1988, Stagg & 
Bailey 1989, Wetherill 1991) have instead adopted various 
aphelion distances Q in the range 8-10 au to distinguish 
Jupiter-family comets from Halley types, while a few authors 
(e.g. Nakamura & Kurahashi 1998) have simply considered 
all comets with P<1000yr as 'periodic'. 

Contrary to the view expressed by Fernandez (1994), 
therefore, a classification based on orbital period alone is 
not very useful. The range of definitions used by different 
investigators to describe essentially the same cometary sub- 
groups appears greatly, and unnecessarily, to complicate any 
comparison of observational or theoretical results obtained 
by different authors. We hope that our historical summary 
has convinced the reader that there is a clear need to improve 
the taxonomy for comet-like bodies in the Solar system! 

3.2 Taxonomy of comets 

Jupiter's zone of control extends between approximately 4 
and 6.6 au. In our classification scheme, a 'comet' has a per- 
ihelion g < 4 au, and so lies formally outside Jupiter's zone 
of control. Nonetheless, Jupiter dominates the lives of all the 
short-period comets. For example, IP/Halley has perihelion 
within the orbit of Venus and aphelion just beyond the or- 
bit of Neptune, yet Jupiter plays the most important role 
in its evolution (Carusi et al. 1988). We propose a four-fold 
division of comets with perihelion distances less than ap- 
proximately 4 au, to include Encke (E), short-period (SP), 
intermediate-period (I) and long-period (L) types. The divi- 
sion is based on aphelion distance, with - as usual - further 
sub-division into the Tisserand parameter classes (Table Q 
denoted by Roman-numeral subscripts. Further embellish- 
ments on this basic scheme could be introduced to accommo- 
date the resonant and non-resonant main-belt asteroids (e.g. 
the Hildas, in the 3:2 mean-motion resonance with Jupiter), 
and objects on Earth-crossing or Earth-approaching orbits 
with periods much less than that of Jupiter, such as near- 
Earth objects. The approximate boundaries between the E, 
SP, I and L comet types are listed in Table |5| 

Working outwards from the Sun, we first come to the 



Object Perihelion 



Aphelion 



E 
SP 
I 
L 



g < 4 
g < 4 
g < 4 
q < 4 



Q<4 
4 < Q < 35 
35 < Q < 1000 
Q > 1000 



J 


4< 


9 


^ 6.6 


Q < 6.6 


JS 


4< 


9 


< 6.6 


6.6<Q<12.0 


JU 


4< 


<? 


^ 6.6 


12.0 < Q < 22.5 


JN 


4< 


<? 


< 6.6 


22.5 <Q< 33.5 


JE 


4 < 


9 


< 6.6 


33.5 < Q < 60.0 


JT 


4< 


9 


< 6.6 


Q > 60.0 



Table 2. The upper panel shows the classification scheme 
for comets. We define Encke- type (E), short-period (SP), 
intermediate-period (I), and long-period (L) comets based on per- 
ihelion and aphelion positions. The lower panel shows the scheme 
for objects whose perihelion is under Jupiter's control. The J class 
describes objects for which both perihelion and aphelion are un- 
der Jupiter's control, the JS class for which the perihelion is un- 
der Jupiter's and the aphelion under Saturn's control, and so 
on. In the final two classes, E stands for Edgeworth-Kuiper belt 
and T for the trans-Neptunian region immediately beyond the 
Edgeworth-Kuiper belt, i.e. the 'trans-EK belt'. 



E-type comets. These have aphelion closer to the Sun than 
approximately 4 au, i.e. Q within the main belt. The proto- 
type is 2P/Encke (Q ~ 4.1 au) which lies on an orbit that 
does not approach closer than a few Hill radii to Jupiter, 
making it stable within the inner Solar system for periods 
of time far greater than that of typical comets. (The orbit 
is not, however, stable for indefinitely long periods, as it is 
subject to long-term secular perturbations which drive the 
perihelion into the Sun on a time-scale of ~100 000 yr; see 
Levison & Duncan 1994, Farinella et al. 1994, Valsecchi et 
al. 1995 for details.) Comet 107P/ Wilson-Harrington, with a 
similar orbit, is a largely inactive body also known as aster- 
oid (4015), and noted by Marsden (1992) to be identical to 
the low-activity comet C/1949 Wl (=1949g) following Bow- 
ell's identification of precovery observations of the asteroid 
on Palomar Schmidt images. Another object, (2201) Oljato, 
has displayed out gassing in the past, and although currently 
classified as an Apollo-type asteroid could well be an evolved 
Encke-type comet (e.g. Weissman et al. 1989, McFadden et 
al. 1993). Many similar bodies in this class appear to be bona 
fide asteroids, rocky bodies perhaps originating via collisions 
in the main asteroid belt. 

Next we come to the SP comets, which include most ob- 
jects conventionally thought of as Halley-type and Jupiter- 
family comets. Many authors have used the Tisserand pa- 
rameter with respect to Jupiter to separate these classes. 
Levison (1996) suggested that comets with Tj < 2.0 should 
be defined as Halley-types and those with Tj > 2.0 should 
be Jupiter-family. However, the precise value of Tj separat- 
ing the two types of behaviour is somewhat arbitrary and 
may depend on period. For example, Rickman & Froeschle 
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(1988) and Bailey (1992) identified a value of the order of 
2.5 as a better marker of a dynamically significant sepa- 
ration. Rather than worry about the precise boundary be- 
tween Jupiter-family and Halley-types, we use the four-fold 
division according to Tisserand parameter introduced in 
Section |T21 namely, Class 1 with Tj < 2.0, Class II with 
2-0 < Tj < 2.5, Class III having 2.5 < T.^ < 2.8 and Class 
IV having Tj > 2.8. Very roughly, these can be thought of 
as corresponding to Halley-types, transition- types, loosely 
bound Jupiter-family comets and tightly bound Jupiter- 
family comets respectively. We call all objects with g < 4 au 
and 4 < Q < 35 au short-period (SP) objects, and differenti- 
ate between objects within the SP class using the Tisserand 
parameter with the sub-class denoted by a Roman numeral 
subscript. So, for example, an SPi object corresponds to a 
Halley-type comet under the conventional classification. 

The third class consists of I-type comets, which have 
perihelion distances less than 4 au and aphelia between ap- 
proximately 35 and 1000 au. Comets with aphelion distances 
less than about 35 au frequently librate around high-order 
resonances with Jupiter (Carusi et al. 1987a, b, 1988; Cham- 
bers 1997) and show evidence for systematic secular changes 
of their orbital elements (Bailey & Emel'yanenko 1996). Re- 
cent examples of I-type comets are C/1995 Ol (Hale-Bopp) 
and 153P/Ikeya-Zhang. They have Tisserand parameters 
with respect to Jupiter of 0.879 and 0.040 respectively, and 
so both fall into the the Ii category. 

The fourth class is the L-type comets, which have per- 
ihelion distances less than 4 au and aphelia greater than 
1000 au. Comets classified under the conventional system as 
long-period fall into either the L or the I class in our system. 
L-types - such as C/1910 Al (the Great January comet) - 
have made no more than their first few passages through 
the inner Solar system. By contrast, I-types are expected to 
have made rather more perihelion passages on average. It is 
sometimes useful to sub-divide L-types into new and young 
(e.g. Fernandez 1981; Wiegert & Tremaine 1999; Horner & 
Evans 2002). New comets, with a > 10 000 au, are first time 
entrants to the Solar system from the Oort Cloud. 

We are of course aware that many comets display ac- 
tivity beyond 4 au from the Sun, and similarly that there 
is an increasing number of known objects with perihelion 
distances beyond 4 au that have been classified as active 
comets (e.g. C/2002V2 [LINEAR] and C/2003A2 [Glea- 
son]) or which have distinctly cometary orbital character- 
istics (e.g. 2002 VQ94). We seek to develop a coherent clas- 
sification scheme that is useful for all comet-like bodies in 
the Solar system. 

From now on, as we move outwards in perihelion, we 
are able to exploit the idea that comet-like bodies lie under 
the control of the giant planets. After the 'comets', the next 
class of objects we encounter are those whose perihelion lies 
within the zone of control of Jupiter. In the conventional 
scheme, these bodies are still classified either as Jupiter- 
family comets, or as asteroids if they show no outgassing. 
In our picture, the designations are listed in Table |5| Here, 
and henceforth, our classification scheme proceeds with the 
first letter designating the planet controlling the perihelion 
and the second letter designating the planet controlling the 
aphelion or the region in which the aphelion lies. The Tis- 
serand parameter is defined with respect to the planet which 
controls the perihelion. So, for instance, a JNiv object has a 



perihelion under Jupiter's control, an aphelion between 22.5 
and 33.5 au and so under Neptune's control, and a Tisserand 
parameter with respect to Jupiter > 2.8; while a JEm type 
has an aphelion between 33.5 and 60 au and so lying in 
the Edgeworth-Kuiper belt and a Tisserand parameter with 
respect to Jupiter in the range 2.5 to 2.8. 

An example of a Jiv object is 29P/Schwassmann- 
Wachmannl, which has a near-circular orbit just beyond 
Jupiter. Its Tisserand parameter is Tj = 3.0 and so it is not 
quite Jupiter-crossing. The object itself is unusual in a cou- 
ple of respects. The low eccentricity of its orbit coupled with 
its unusually great size (~20km radius) allows it to be fol- 
lowed with ease through its entire orbit. It undergoes regular 
outbursts, with its visual magnitude rising from its usual 17 
(at perihelion) by a number of magnitudes, reaching as high 
as 10 at its brightest. 

An example of a JSiv type object is 39P/Oterma. It 
was originally discovered in 1942 orbiting in a Hilda-type 3:2 
mean-motion resonance with Jupiter, with a low-eccentricity 
orbit and a period of a mere 8 years. At that time, its orbit 
lay between those of Mars and Jupiter and it would have 
been classified as an SPiv comet in our picture (although 
clearly as a resonant 'Hilda' in a more complete scheme). 
However, it was later noted that it had entered the 3:2 res- 
onance in 1937 from the 2:3 resonance, and returned to the 
2:3 in 1963 (e.g., Kazimirchak-Polanskaya 1972). Continuous 
periods in any one resonance are short owing to the orbit's 
unstable nature, but it is noteworthy that Tj ~ 3.0 was 
approximately conserved during this evolution. The present 
orbit lies outside that of Jupiter, with a current perihelion 
distance of 5.46 au and an aphelion of 9.02 au. This object 
seems to be a good example of a recent 'handing down' of a 
Centaur from Saturn's control to Jupiter, with further peri- 
ods as an SP comet seeming quite probable. Finally, recent 
discoveries like C/2002 Al (LINEAR) and C/2002 A2 (LIN- 
EAR), with semi-major axes close to that of Uranus, furnish 
examples of JNni objects. 



3.3 Resonances 

A key feature of the orbital evolution of Centaurs and many 
other comet-like objects in the outer Solar system is the 
tendency for the semi-major axis to lie close to one or an- 
other mean- motion resonance with a planet. This behaviour 
is exemplified by the so-called Plutinos (objects in the 2:3 
mean-motion resonance with Neptune, i.e. around a ~ 39.5 
au), which appear clearly in Fig. |5| and also by the evolu- 
tion of 2000 FZ53 shown in the top left plot of Fig. The 
occurrence or otherwise of a long-lived resonance will clearly 
affect the transfer probability from one part of the (q, Q, Tp) 
phase space to another and is an important consideration 
in the dynamical evolution of outer Solar system objects. 
There are recently discovered objects with semi-major axes 
close to essentially all the significant Neptunian resonances, 
from 1:1 [e.g. 2001 QR322 and 2002 PN34] to 1:8 [e.g. (54520) 
2000 PJ30]. However, detailed discussion of the classification 
and role of such resonances is beyond the scope of this paper. 
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Figure 2. Plot of eccentricity versus semi-major axis showing the boundaries of the classification scheme. The red circles are listed as 
Centaurs and Scattered Disk Objects by the Minor Planet Center, while the green circles are Edgeworth-Kuiper belt members. In both 
cases, we only record objects if their observed arcs are 30 days or more. 
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Figure 3. The evolution of dynamical quantities for a clone of the Centaur 2000 FZ53. The left panel of plots show the variation of 
semi-major axis, perihelion, aphelion, eccentricity, inclination and Tisserand parameter (in this case with respect to Uranus) against 
time in years. The rightmost plot shows the evolution in the plane of semi-major axis and eccentricity, with the different categories of 
the classification scheme marked. The clone spends most of the first 1 Myr as a UE object and ends in the NE category. Notice the very 
sharp transition from UN to UE object at ~ 2 Myr. [Starting orbital elements of the clone are a = 23.595 au, e = 0.497, i = 34.865° at 
JD 2451120.5.] 



4 CENTAURS AND TRANS-NEPTUNIAN 
OBJECTS 



Table 121 gives the classification scheme for objects under the 
control of giant planets other than Jupiter. As before, the 
notation is that the first letter designates the planet control- 
ling the perihelion, and the second letter the planet control- 
ling the aphelion or the region in which the aphelion lies. 
So, a UN object has a perihelion controlled by Uranus and 
an aphelion controlled by Neptune, and so on. The letters 
EK denote objects with perihelia and aphelia lying close 
to or within the classical Kuiper or Edgeworth-Kuiper belt. 
Trans-EK belt or T objects, on the other hand, may have 
much larger aphelion distances, and (perhaps depending on 
their perihelion distances) may be stable 'outer' objects as 
described by Emel'yanenko et al. (2003) , or dynamically ac- 
tive objects similar to those conventionally known as scat- 
tered disc objects. Both EK and T classes have perihelia 
Table 3. Object classification for the trans-Jovian region. The beyond Neptune's control (33.5 au) but less than the con- 

first letter designates the planet controlling the perihelion, the ventional outer edge of the Edgeworth-Kuiper beh (~ 60 

second letter the planet controlling the aphelion or the region au). They differ in that the aphelion Q is less than approx- 

in which the aphelion lies, with the final two classes EK and T imately 60 au for the former and greater than 60 au for the 

being beyond all the giant planets. (S = Saturn, U = Uranus, N latter. As for the comets, the information on the zone of 

= Neptune, EK = Edgeworth-Kuiper belt, T = trans-EK belt). control must be supplemented with the Tisserand parame- 

ter with respect to the controlling planet to give a complete 
classification. The Tisserand parameter class is given as a 
subscript. It is interesting to note the recent discovery of sev- 
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Figure 4. The evolution of dynamical quantities for a clone of Chiron. The Tisserand parameter plotted here is that with respect to 
Jupiter. The clone starts as an SU object, evolves into a short-period Jupiter-family comet and then is finally ejected from the Solar 
system. Notice how evolution often proceeds along lines of almost constant aphelion or perihelion in the rightmost panel. [Starting orbital 
elements of the clone are a = 13.599 au, e = 0.384, i = 6.879° at JD 2451220.5.] 
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Figure 5. The evolution of dynamical quantities for a clone of 1998 TF35. The Tisserand parameter plotted here is that with respect 
to Jupiter. The clone starts as a UE object, evolves into a short-period Jupiter-family comet and then is ejected from the Solar system. 
[Starting orbital elements of the clone are a = 26.412, e = 0.369, i = 12.564° at JD 2451220.5.] 



eral high-inclination Centaurs, for example 2002 VQ94 and 
2002 XU93, which we classify as STi and UTi respectively. 

The first Centaur discovered, (2060) Chiron, has a per- 
ihelion distance of 8.4 au and an aphelion of 18.8 au. It 
is classified as an SUiv object, as its Tisserand parameter 
Ts = 2.9. This is one of the largest Centaurs and is active at 
perihelion. Hahn & Bailey (1990) showed that Chiron's orbit 
is rapidly evolving and that it may have been a short-period 
comet in the past, and will probably become one again in the 
future. The second Centaur discovered, (5145) Pholus, has 
a perihelion of 8.7 au and an aphelion of 32.1 au, together 
with a Tisserand parameter Ts of 2.6. It is classified as an 
SNiii object. Pholus's orbit is also unstable with a charac- 
teristic lifetime of ^ 10^ yr. As another example, (10199) 
Chariklo has a perihelion of 13.1 au which is just within 
Uranus's control and an aphelion of 18.6 au. Its Tisserand 
parameter with respect to Uranus is 2.9, making its clas- 
sification Uiv. Chariklo's orbit is more stable than either 
Chiron's or Pholus's, as the dynamical lifetime is roughly 
a factor of 10 times longer. Not all orbits in this part of 
the Solar system are evolving. For example, Holman (1997) 
found that a few low-eccentricity objects between 24 and 
27 au (^^0.3% of an initial low-eccentricity primordial sam- 
ple) may survive for 4.5 Gyr without significant dynamical 
evolution. These objects lie in the Niv class. Others (e.g. 
2OOIQR322, also classified as Niv) are recognized as long- 
lived Neptune- Trojans (Marsden 2003), and it is probably 
only a matter of time before similar examples of Trojans are 
confirmed around Uranus as well. 



Almost all the known Centaurs lie within the Tisserand 
parameter classes III and IV. There are only a few exceptions 
- for example, 2000 FZ53 has a Tisserand parameter with re- 
spect to its controlling planet (Uranus) of 2.4. However, the 
nearly isotropic influx from the Oort Cloud must inevitably 
produce at least a few high-inclination Centaurs. Hence, the 
other classes are needed, and their apparent emptiness at 
present is largely a selection effect. 

Fig-IHshows the plane of semi-major axis and eccentric- 
ity. The red circles show the current location of the 52 ob- 
jects listed as Centaurs and Scattered Disk Objects in 2002 
June by the Minor Planet Center'^. A further 16 objects lie 
at semi-major axes too great to be seen on the plot. These 
objects were chosen from the full list of Centaurs by requir- 
ing the observable arc to be greater than 30 days. Objects 
with shorter arcs are excluded as their orbits often change 
substantially from the original estimates when the arc is ex- 
tended. The known population depends both on the true 
numbers in the various regions and on discovery selection 
effects. Orbits under the control of Uranus and Neptune are 
expected to be more long-lived than those under the control 
of Saturn and Jupiter. However, this is counterbalanced by 
the fact that objects nearer the Sun are more easily discov- 
ered. The question can be studied further by modelling the 
efficiency of observational surveys (Jedicke & Herron 1997). 
For the present Centaur population, there is a suggestion in 

^ http:/ /cfa-www. harvard.edu/iau/lists/Centaurs. html 
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Figure 6. The evolution of comets P/1997 T3 (Lagerkvist-Carsenty), upper panels, and P/1998 U3 (Jager), lower panels. On the left, 
the numerically integrated orbit is shown together with the orbits of Jupiter, Saturn and Uranus. On the right, the evolution of the 
comet in the plane of semi-major axis and eccentricity is presented. We have suppressed some of the labels to reduce notational clutter. 
Notice that P/1997 T3 shows a perihelion-aphelion exchange on transference from the Saturn dominated regime to Jupiter family. For 
comet Jager, the evolution proceeds through the encounter at almost constant aphelion. 



Fig. |21of a substantial population in the UE category. These 
may be being scattered down from planet to planet having 
originated in the Edgeworth-Kuiper belt. 

The green circles in Fig.|5|show the current classification 
of those 373 objects with arcs greater than 30 days belonging 
to the Edgeworth-Kuiper belt. Again, our data were down- 
loaded in June 2002 from the Minor Planet Center^. The two 
most populous categories are NE and EK. The Plutinos, a 
population of objects trapped in the 2:3 mean-motion reso- 
nance with Neptune (e.g. Jewitt & Luu 1996), fall into either 
the NE or EK category depending on eccentricity. The reso- 
nance acts to prevent close approaches to Neptune, although 
it is unlikely that such a mechanism would generally protect 
all the objects from eventual encounters with Neptune over 
time-scales comparable to the age of the Solar system. The 
EK category contains most of what are conventionally called 
the 'classical' Kuiper or Edgeworth-Kuiper belt population, 
while the T category contains most of what are convention- 
ally called scattered disc objects. These last two classes can 
be further subdivided (Emel'yanenko et al. 2003). 



5 NUMERICAL EXAMPLES 

We now illustrate the usefulness of the classification scheme 
with some examples drawn from a suite of numerical sim- 
ulations. Clones of the current Centaur population are in- 
tegrated in the forward and backward direction for up to 

6 Myr using the hybrid integrator in the Mercury pro- 
gram (Chambers 1999). The orbits are integrated under the 
gravitational effects of the Sun and the four major planets 
only. The time-step is 120 days. Once objects passed beyond 
1000 au, their orbits are no longer followed. 



5.1 From UEn to NEn object 

Fig-Elshows the evolution of a clone of 2000 FZ53. For most 
of the first 1 Myr its perihelion is under the control of Uranus 
and its aphelion lies in the Edgeworth-Kuiper belt. The Tis- 
serand parameter with respect to Uranus is T\j — 2.4. The 

^ http:/ /cfa-www. harvard.edu/iau/lists/TNOs. html 



clone starts off as an SEn object and moves rapidly into the 
UEii zone. 

During the course of the simulation, the clone is driven 
gradually towards a stable orbit with perihelion close to Nep- 
tune (an NE object). This illustrates one route by which 
an object can approach the Edgeworth-Kuiper belt from 
the Centaur region. There are two sharp transitions visi- 
ble in Fig. 01 The first change, at 2.0 Myr from UN to UE, 
is caused by a close encounter at aphelion with Neptune; 
the second, at 2.9 Myr within the UE category, is caused by 
an encounter at perihelion with Uranus. At the end of the 
6 Myr simulation, the object has a perihelion just beyond 
the orbit of Neptune, a high inclination, and a low eccen- 
tricity. The value of its Tisserand parameter with respect 
to Neptune is Tn — 2.4. The orbit is similar to those of a 
number of objects currently residing in the inner regions of 
the Edgeworth-Kuiper belt, such as 1999CPi33. 

The stability of the semi-major axis towards the end 
of the 2000 FZ53 simulation suggests that it lies in a mean- 
motion resonance with Neptune. This is borne out by cal- 
culation of the orbital period of the clone, which is ~202 yr, 
suggesting that the object probably lies in the 4:5 mean- 
motion resonance with Neptune (period ~164yr). It is pos- 
sible that the orbit will continue to drift outward and sta- 
bilize after the end of the simulation, illustrating a rather 
rare example of outward rather than inward evolution. The 
ultimate destiny of such an object may therefore be to join 
the Edgeworth-Kuiper belt or the region beyond it. 

5.2 From SUiv object to SPiv comet 

Fig. 21 shows the evolution of a clone of Chiron (cf. Hahn & 
Bailey 1990). This object illustrates one of the supply routes 
for the short-period comets. The clone starts with its peri- 
helion under the control of Saturn, but its aphelion wanders 
between the control of Saturn, Uranus and Neptune respec- 
tively for the first 2x10* yr. After this, the clone falls under 
Jupiter's control through a succession of aphelion-perihelion 
interchanges. The Tisserand parameter Tj falls below 3 and 
the clone spends the next part of its lifetime (~ 7x 10* yr) 
as a short-period comet before leaving the SP region. This 
is longer than the typical fading time (~ 1.2 x 10* yr) for 
Jupiter-family comets found by Levison & Duncan (1997). 
It is noteworthy that this 200 km diameter object has a sig- 



Populations of Comet-Like Bodies 9 



nificant probability of evolving to an orbit with perihelion 
distance less than that of the Earth (Hahn & Bailey 1990). 
This emphasizes the potential significance for the inner so- 
lar system of so-called 'giant' comets and their disintegration 
products when time-scales of the order of 10^ yr or more are 
considered (e.g., Napier 2001). From the characteristic shape 
of the trajectory in the rightmost panel of Fig. 21 we see that 
the evolution proceeds along lines of roughly constant per- 
ihelion or aphelion both before and after the object comes 
under Jupiter's control. The part of the clone's life when it 
is controlled by Jupiter is characterized by it being highly 
chaotic and showing rapid dynamical evolution. 

Note that the clone starts with a Tisserand parameter 
with respect to Saturn Ts — 2.9. During its lifetime as a 
short-period comet, its Tisserand parameter with respect 
to Jupiter Tj = 2.8. Just as in our previous example, the 
Tisserand parameter with respect to the controlling planet 
is nearly preserved during the whole evolution. 

The fact that the clone stays in an area in which it 
will be active as a comet for such a long time suggests that 
cometary bodies can be captured into the inner Solar system 
for long enough to have a reasonable chance of decoupling 
from Jupiter, whether by non-gravitational forces such as 
outgassing (cf. Kresak 1982, Harris & Bailey 1998, Asher, 
Bailey & Steel 2001) or collisional mechanisms. Integrations 
occasionally provide examples of Encke-type orbits evolving 
to become Jupiter-family short-period comets, but the re- 
verse evolution by gravitational means is so rare that it has 
never been seen (Carusi & Valsecchi 1987). It is interesting 
to note that the stability of an orbit approaching Saturn 
is lower than one approaching merely Uranus or Neptune. 
Similarly, an orbit in the vicinity of Jupiter is yet more un- 
stable. This adds weight to our argument that classifying all 
objects in the trans- Jovian region merely as 'Centaurs' is in- 
sufficient to explain the full richness of dynamical behaviour 
in the region. 



5.3 From UEiv object to SPiv comet 

Fig.|S]shows the evolution of a clone of 1998 TF35. It starts 
off as a UE object with Tisserand parameter Ttj — 2.9. 
This object is of interest since it passes through a num- 
ber of different categories (UE, UN, NE, N, U, SU and S) 
before finally becoming a short-period comet and then suf- 
fering ejection from the Solar system. Generally, it moves 
through the categories sequentially rather than frequently 
jumping back and forth between classes. The clone begins 
in an orbit influenced both by Neptune and Uranus, and 
distant encounters provide gradual perturbations for most 
of its lifetime. Eventually, encounters with Neptune move 
the object into a near-circular orbit controlled only by that 
planet. This period of evolution ends with a close encounter 
which moves the clone's perihelion back to the control of 
Uranus. The perihelion drops further until the object en- 
counters Saturn, at which point it is rapidly handed down 
to Jupiter's control, where it becomes a short-period comet 
for a brief time (~ 10* yr) before being ejected from the So- 
lar system. Whilst a short-period comet, Tj = 2.9, so there 
is again excellent conservation of the Tisserand parameter 
with respect to the controlling planet. 



5.4 Comets Lagerkvist-Carsenty and Jager 

As a final example, let us consider two objects convention- 
ally classified as short-period Jupiter-family comets, namely 
P/1997T3 (Lagerkvist-Carsenty) and P/1998U3 (Jager). 
Numerical integrations by Lagerkvist et al. (2000) showed 
that both objects suffered recent close encounters with Sat- 
urn that drastically changed their orbits. Fig. |S] shows the 
evolution of the orbits of both comets over the last century. 
Comet Lagerkvist-Carsenty started out as an SUiv object 
in January 1900. A close encounter with Saturn in October 
1954 transferred its osculating elements through almost in- 
stantaneous SN and SE phases into its present JSiv orbit. 
This proceeded through a perihelion-aphelion interchange, 
with the old perihelion becoming the new aphelion. By con- 
trast. Comet Jager started as an Sm object in January 1900 
and underwent a direct transition to its present SPu orbit in 
July 1991. The evolution proceeded at nearly constant aphe- 
lion, as can be deduced both from the plot of the orbit and 
from its evolutionary track in the plane of semi-major axis 
and eccentricity. In this case, the close approach caused the 
Tisserand parameter class to change, the controlling planet 
switching from Saturn to Jupiter. 

These orbital integrations provide illustrations of why 
the classification scheme is useful. The scheme is based 
around the aphelion and perihelion distances. Close encoun- 
ters often preserve aphelion or perihelion to lowest order, or 
a perihelion-aphelion interchange may take place. Very of- 
ten, the Tisserand parameter with respect to the controlling 
planet is conserved to a good approximation as well. 



6 CONCLUSIONS 

The main aim of the paper is to introduce a new classifi- 
cation system for comet-like bodies. Minor bodies between 
Saturn and Neptune are often described simply as 'Centaurs' 
and those beyond Neptune simply as 'Kuiper Belt Objects'. 
This is not very enlightening as the histories and fates of 
such bodies may be very different. For example, simulations 
show that some Centaurs are in long-lived and stable orbits, 
whereas others are dynamically highly active and will evolve 
rapidly. 

The main focus of this paper has been on the Centaurs 
and short-period comets. Our proposition is to classify the 
comet-like objects beyond Jupiter according to the planets 
that control the evolution of their perihelion and aphelion. 
For example, an SN object has a perihelion under the con- 
trol of Saturn and an aphelion under the control of Nep- 
tune, while a UE object has a perihelion under the control 
of Uranus and an aphelion in the Edgeworth-Kuiper belt. 
This provides 20 dynamically distinct categories of outer So- 
lar system objects in the Jovian and trans- Jovian regions. 
The evolutionary tracks of bodies like Centaurs often show 
periods in which the aphelion or perihelion distances are 
individually rather well conserved, or encounters in which 
the old perihelion becomes the new aphelion or vice versa 
('perihelion-aphelion interchanges') . 

The objects with smaller perihelion distances, which 
we designate as 'comets', are dominated by Jupiter and 
not by the planets nearest at perihelion or aphelion. In our 
scheme, comets have perihelion distance less than 4 au. They 
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are subdivided into Encke-type, short-period, intermediate- 
period and long-period according to aphelion distance. Fol- 
lowing a succession of authors beginning with Kresak (1972), 
we favour further sub-divisions based on the Tisserand pa- 
rameter. Specifically, we subdivide the categories of comets 
into: Class I which has Tj < 2.0, Class II which has 
2-0 < Tj < 2.5, Class III having 2.5 < Tj < 2.8 and Class 
IV having Tj > 2.8. Very roughly, these can be thought of 
as corresponding to Halley-types, transition-types, loosely 
bound Jupiter-family comets and tightly bound Jupiter- 
family comets respectively. This idea is then extended to all 
comet-like bodies by identifying the Tisserand parameter of 
the planet controlling the perihelion with a subscript. 

Given the aphelion and perihelion distance, together 
with the Tisserand parameter of the planet controlling the 
perihelion, it is straightforward to find the instantaneous 
classification of any Solar system object. Our new classifi- 
cation scheme extends the existing taxonomy for comets to 
cover all comet-like bodies in the Solar system. 



ACKNOWLEDGMENTS 

This research was supported by the Particle Physics and As- 
tronomy Research Council (JH), the Royal Society (NWE) 
and the Northern Ireland Department of Culture, Arts and 
Leisure (MEB, DJA). DJA acknowledges the hospitality of 
the Japan Spaceguard Association during work on this pa- 
per. We thank the referee for very constructive comments. 



REFERENCES 

Applcgate J.H., Douglas M.R., Giirsel Y., Sussman G.J., Wisdom, 

J., 1986, AJ, 92, 176 
Arnold V.I., Kozlov V.V., Neishtadt A.I., 1987 in Arnold V.I., 

ed.. Dynamical Systems III. Springer Verlag, New York, p. 70 
Asher D.J., Bailey M.E., Hahn G., Steel D.I., 1994, MNRAS, 267, 

26 

Asher D.J., Bailey M.E., Steel D.I., 2001, in Marov M., Rickman 
H., eds, Astrophys. Space Sci. Libr. 261, CoUisional Processes 
in the Solar System. Kluwer, Dordrecht, p. 121 
Bailey M.E., 1992, Cel. Mech. Dyn. Astron., 54, 49 
Bailey M.E., Emel'yanenko V.V., 1996, MNRAS, 278, 1087 
Carusi A., Valsecchi G.B., 1987, European Regional Astronomy 

Meeting of the lAU, Volume 2, 21 
Carusi A., Kresak L., Perozzi E., Valsecchi G.B., 1987a, European 

Regional Astronomy Meeting of the lAU, Volume 2, 29 
Carusi A., Kresak L., Perozzi E., Valsecchi G.B., 1987b, A&A, 
187, 899 

Carusi A., Valsecchi G.B., Kresak L., Perozzi E., 1988, Cel. Mech., 
43, 319 

Chambers J.E., 1997, Icarus, 125, 32 

Chambers J.E., 1999, MNRAS, 304, 793 

Charnoz S., Thebault P., Brahic A., 2001, A&A, 373, 683 

Crovisier J., 2001, in Murdin P., ed.. Encyclopedia of Astronomy 

and Astrophysics. lOP Publishing Ltd and Nature Publishing 

Group, p. 446 

Dones L., Levison H.F., Duncan M., 1996, in Rettig T.W., Hahn 
J.M., eds, ASP Conf. Ser. 107, Completing the Inventory of 
the Solar System. ASP, San Francisco, p. 233 
Duncan M.J., Levison H.F., 1997, Sci, 276, 1670 
Duncan M., Quinn T., Tremaine S., 1988, ApJ, 328, L69 
Emel'yanenko V.V., Asher D.J., Bailey M.E., 2003, MNRAS, 338, 
443 



Everhart E., 1972, Astrophysical Letters, 10, 131 
Everhart E., 1973, AJ, 78, 329 

Everhart E., 1974, in Cristescu C, Klepczynski W.J., Milet B., 
eds, Proc. lAU Coll. 22, Asteroids, comets, meteoric matter, 
p. 223 

Everhart E., 1977, in Delsemme A.H., ed., Proc. lAU Coll. 39, 
Comets Asteroids Meteorites: Interrelations, Evolution and 
Origins. University of Toledo, Toledo, p. 99 

Farinella P., Proeschle Ch., Froeschle C, Gonczi R. , Hahn C, 
Morbidelli A., Valsecchi G.B., 1994, Nat, 371, 314 

Fernandez J. A., 1980 MNRAS, 192, 481 

Fernandez J.A., 1981, A&A, 96, 26 

Fernandez J. A., 1985, Icarus, 64, 308 

Fernandez J. A., 1994, in Milani A., Di Martino M., Cellino A., 
eds, Proc. lAU Symp. 160, Asteroids, Comets, Meteors 1993. 
Kluwer, Dordrecht, p. 223 

Fernandez J.A., Gallardo T., 1994, A&A, 281, 911 

Fernandez J.A., Gallardo T., Brunini A., 2002, Icarus, 159, 358 

Green D.W.E., 2000, International Comet Quarterly, 22, 2 

Hahn C, Bailey M.E., 1990, Nat, 348, 132 

Harris N.W., Bailey M.E., 1998, MNRAS, 297, 1227 

Holman M.J., 1997, Nat, 387, 785 

Horner J., Evans N.W., 2002, MNRAS, 335, 641 

Jcdicke R., Herron J.D., 1997, Icarus, 127, 494 

Jewitt D., Luu J.X., 1993, Nat, 362, 730 

Jewitt D., Luu J.X., 1996, in Rettig T.W., Hahn J.M., eds, ASP 
Conf. Ser. 107, Completing the Inventory of the Solar System. 
ASP, San Francisco, p. 255 

Kazimirchak-Polonskaya E.I., 1972, in Chebotarev G.A., 
Kazimirchak-Polonskaya E.I., Marsden B.C., eds, Proc. lAU 
Symp. 45, The Motion, Evolution of Orbits, and Origin of 
Comets. Reidel, Dordrecht, p. 373 

Kowal C.T., 1979, in Gehrels T., ed.. Asteroids. University of 
Arizona Press, Tucson, p. 436 

Kresak L., 1972, in Chebotarev G.A., Kazimirchak-Polonskaya 
E.I., Marsden B.C., eds, Proc. lAU Symp. 45, The Motion, 
Evolution of Orbits, and Origin of Comets. Reidel, Dordrecht, 
p. 503 

Kresak L., 1980, The Moon and Planets, 22, 83 

Kresak L., 1982, in Fricke W., Teleki C, eds. Sun and Planetary 

System. Reidel, Dordrecht, p. 361 
Kresak L., 1983, in West R.M., ed.. Highlights of Astronomy, 6, 

377 

Kresak L., 1985, in Carusi A., Valsecchi G.B., eds, Proc. lAU 
Coll. 83, Dynamics of Comets: Their Origin and Evolution. 
Reidel, Dordrecht, p. 279 

Larsen J. A. et al., 2001, AJ, 121, 562 

Lagerkvist C.-L, Hahn G., Karlsson O., Carsenty U., 2000, A&A, 
362, 406 

Levison H.F., 1996, in Rettig T.W., Hahn J.M., eds, ASP Conf 
Ser. 107, Completing the Inventory of the Solar System. ASP, 
San Francisco, p. 173 

Levison H.F., Duncan M.J., 1994, Icarus, 108, 18 

Levison H.F., Duncan M.J., 1997, Icarus, 127, 13 

Luu J.X., Jewitt D., 1990, AJ, 100, 913 

Luu J.X., Jewitt D., Trujillo C.A., Hergenrother C.W., Chen J., 

Offutt W.B., 1997, Nat, 387, 573 
Manara A., Valsecchi G.B., 1991, A&A, 249, 269 
Marsden B.C., 1992, lAU Circ. 5585 

Marsden B.C., 1999a, http:/ /cfa- www. harvard.edu/cfa/ps/pressinfo/200TNOs. html 

Marsden B.C., 1999b, Minor Planet Circ. 35155 

Marsden B.C., 2003, Minor Planet Elect. Circ. 2003- A55 

Marsden B.C., Williams, G.V., 1999, Catalogue of Cometary Or- 
bits 1999. Minor Planet Center, Cambridge, Massachusetts 

McFadden L.A., Cochran A.L., Barker E.S., Cruikshank D.P., 
Hartmann W.K., 1993, J. Geophys. Res., 98, 3031 

Murray CD., Dermott S.F., 1999, Solar System Dynamics. Cam- 
bridge University Press, Cambridge 



Populations of Comet-Like Bodies 11 

Nakamura T., Kurahashi, H., 1998, AJ, 115, 848 

Napier W.M. 2001, MNRAS, 321, 463 

Quinn T., Tremaine S., Duncan M., 1990, ApJ, 355, 667 

Rickman H., Proeschle CI., 1988, Cel. Mech., 43, 243 

Scotti J.V., 1992, lAU Circ, 5434 

Stagg C.R., Bailey M.E., 1989, MNRAS, 241, 507 

Stern A., Campins, H., 1996, Nat, 382, 507 

Vaghi S., 1973a, A&A, 24, 107 

Vaghi S., 1973b, A&A, 29, 85 

Valsecchi G.B., Morbidelli A., Gonczi R., FarincUa P., Frocschle 
Ch., Proeschle C, 1995, Icarus, 118, 169 

Weissman P.R., 2001, in Murdin P., ed.. Encyclopedia of Astron- 
omy and Astrophysics, lOP Publishing Ltd and Nature Pub- 
lishing Group, p. 1930 

Weissman P.R., A'Heajrn M.F., McFadden L.A., Rickman H., 
1989, in Binzel R.P., Gehrels T., Matthews M.S., eds. As- 
teroids II, University of Arizona Press, Tucson, p. 880 

Wcthcrill G.W., 1991, in Ncwburn R.L. Jr, Ncugebaucr M., Rahe 
J., cds, Proc. lAU Coll. 116, Comets in the Post-Halley Era, 
Vol. 1. Kluwcr, Dordrecht, p. 537 

Whipple F.L., 1978, The Moon and Planets, 18, 343 

Wiegert P., Tremaine S., 1999, Icarus, 137, 84 



This figure "Chiron.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0304319vl 



This figure "Comets.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0304319vl 



This figure "Currentobs.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0304319vl 



This figure "TF35.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0304319vl 



This figure "FZ53.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0304319vl 



